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N-flation is a promising embedding of inflation in string theory in which many string 
axions combine to drive inflation [||. We characterize the dynamics of a general N-flation 
model with non-degenerate axion masses. Although the precise mass of a single axion 
depends on compactification details in a complicated way, the distribution of masses can 
be computed with very limited knowledge of microscopies: the shape of the mass distri¬ 
bution is an emergent property. We use random matrix theory to show that a typical 
N-flation model has a spectrum of masses distributed according to the Marcenko-Pastur 
law. This distribution depends on a single parameter, the number of axions divided by the 
dimension of the moduli space. We use this result to describe the inflationary dynamics 
and phenomenology of a general N-flation model. We produce an ensemble of models and 
use numerical integration to track the axions’ evolution and the resulting scalar power 
spectrum. For realistic initial conditions, the power spectrum is considerably more red 
than in single-field m 1 2 (j) 2 inflation. We conclude that random matrix models of N-flation 
are surprisingly tractable and have a rich phenomenology that differs in testable ways from 
that of single-field m 2 (j) 2 inflation. 
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1. Introduction 

Inflationary model-building in string theory [|2|,[3j is advancing rapidly, impelled by the 
strong and growing evidence for inflation, and facilitated by a series of technical advances 
|J,[5],|Sj in string theory. Among the inflationary models that have remained difficult to 
embed in string theory are those known as large-held models, in which the inhaton vacuum 
expectation value (vev) varies by much more than one Planck mass during inflation j?|. 
This class includes what are arguably the simplest inhationary potentials, those quadratic 
IP or quartic in the inhaton. 

The controllability of these models is a subject of debate 0, and depends strongly 
on the relation between particle physics and gravity near the Planck scale. Many authors 
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take the view that configurations with trans-Planckian vevs are controllable provided that 
the energy density is well below the Planck scale. Because large-field inflation requires 
a very weakly coupled field, potentials such as Xcj/ 1 and m 2 (j) 2 automatically have safely 
sub-Planckian energy densities. This argument is reliable if we restrict our attention to 
quantum corrections that involve higher-order terms in the spacetime curvature 1Z , since 
TZ -C M 2 during the portion of inflation that is observationally relevant. In a model that 
does not unify gravity with the other forces we might ignore gravitationally induced correc¬ 
tions to the particle physics effective potential, and only worry about quantum corrections 
to the spacetime background. In this context, trans-Planckian vevs can indeed be benign. 

However, gravitationally-suppressed couplings between the fields that contribute to 
T //(/ generically appear in models that unify gravity and particle physics, such as super¬ 
gravity and superstring models. The problem is that loop corrections generically become 
large when vevs, instead of energy densities, exceed the Planck scale; for example, it is very 
difficult to forbid corrections to the Kahler potential that are controlled by field vevs in 
Planck units. Experience with supergravity and superstring models strongly suggests that 
Planck scale vevs usually imply large corrections to the inflaton effective potential. This 
observation has been used to put bounds on inflationary observables, such as the tensor 
power spectrum 0].l 

Our view is that because trans-Planckian vevs are always dangerous, if not necessar¬ 
ily fatal, in string theory and supergravity models, it is worthwhile to construct specific 
models in which these problems can be overcome. In the best case, understanding trans- 
Planckian vevs in a concrete string theory construction would give a useful toy model for 
understanding general quantum gravity corrections to large-field models. 

Dimopoulos, Kachru, McGreevy, and Wacker [jl| have proposed a clever solution to this 
set of problems.i Combining IV 3> 1 fields with sub-Planckian displacements into a single 
effective inflaton, they obtained m 2 (j) 2 inflation without having to build a trans-Planckian 
single-field displacement. The resulting model, which they called N-flation, is an example of 
assisted [ p~5| , p~3| , pL6| , pr7| , p~8|l inflation. They further explained that N-flation can be embedded 
in string theory by choosing the N fields to be axions. Axions are typically present in 
large numbers in string compactifications, and even when all other moduli are stabilized, 


4 There are a few exceptions to this rule [11,12]. 

4 For earlier work in this direction, see [13]. For an application of these ideas to quintessence, 


see ]|14j]. 


2 








the axion potentials remain rather flat as a consequence of well-known nonrenormalization 
theorems m- 

The axions are pseudo-Nambu-Goldstone bosons associated with nonperturbatively- 
broken shift symmetries. Hence, N-flation may be categorized as an approximately shift- 
symmetric, closed string, large-field model.I Shift-symmetric D-brane configurations have 
been suggested [Q as a natural origin of flat inflaton potentials in string theory; however, 
the effects that break such open string shift symmetries |22[] are often too strong to leave a 
suitably flat potential. (See also [^3| for a general explanation of the role of approximate 
shift symmetries in preserving flat potentials.) 

When the individual axion displacements are small in Planck units, the leading term 
in the potential is a mass matrix that connects the N axion fields. Choosing a basis 
in which the kinetic terms are canonical and the mass terms are diagonal, the data of 
an N-flation model becomes simply a list of N masses. The authors of J|!] focused their 
attention on the special case in which the masses are identical, and showed that in that 
case the observational predictions are precisely those of m 2 © 2 chaotic inflation. They 
argued that inflation still occurs in the general case with non-degenerate masses, with 
qualitatively similar results. The dynamics is somewhat more complicated, however, as 
the most massive axions relax first to their minima, diminishing the Hubble friction and 
causing the lighter axions to roll more rapidly. 

The purpose of this note is to give a quantitative treatment of this general case of 
N-flation, in which all the masses are distinct. Two natural questions arise in this context. 
First, what are the dynamics and the observational predictions of an N-flation model with a 
given set of masses rrq? In particular, how does a spread in masses affect the spectral index 
n s ? Second, what spectrum of masses should one expect from microscopic considerations? 
The mass of each axion depends on a formidable array of compactification details. In the 
context of KKLT moduli stabilization ||, these details include some as-yet uncomputed 
factors, such as fluctuation determinants for Euclidean D3-branes, in addition to better- 
known but highly unwieldy expressions involving the fluxes and the Kahler potential. 
Computing the precise masses of N axions in a nontrivial, stabilized compactification 
would therefore appear to be extremely difficult. Might we nevertheless be able to deduce 
something about the typical spacing of axion masses? 


Natural inflation 120,211, in contrast, is a small-field model driven by a single axion. 
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Surprisingly, we find an essentially universal result for the axion masses in N-flation: 
the probability distribution for the masses-squared is given by an analytic formula known 
as the Marcenko-Pastur law. The shape of this curve depends on a single parameter (3 
determined by the dimensions of the Kahler and complex structure moduli spaces. We 
describe this result as universal because it does not depend on specific details of the 
compactification, such as the intersection numbers, the choice of fluxes, or the location in 
moduli space. The result is also insensitive to superpotential corrections such as threshold 
corrections and instanton determinants. The shape of the mass distribution depends only 
on the basic structure of the mass matrix, which is specified by the supergravity potential. 
(The overall mass scale is not determined a priori ; as in [[IJ, we fix it by requiring the 
observed amplitude of density perturbations.) 

This welcome situation arises because of the size (typically larger than 200 x 200) of 
the matrices involved: the distribution of the eigenvalues of random matrices is most easily 
characterized when the matrices are large. In this way, the simplifications enjoyed by large 
random matrices [p4| , |25|j2^ , p7| help N-flation to generate a prediction for the scalar spec¬ 
tral index which is substantially (though not completely) independent of compactification 
details. This is a striking example of a problem that is tractable because of, rather than 
in spite of, the high dimensionality of the scalar field space in string theory. 

Equipped with the Marcenko-Pastur law for the mass distribution, we numerically 
integrate the equations of motion for a large sample of realizations of N-flation. We find 
that for a range of initial conditions and for most values of f3, the scalar power spectrum is 
typically more red (i.e. farther from scale-invariant) than in the case where all the masses 
are identical. 

In addition to developing a random matrix theory calculation of the axion mass spec¬ 
trum, we treat the inflationary dynamics and perturbation spectrum produced by N-flation. 
(For related earlier work, see [[28|. |20l| . Additionally, Kaloper and Liddle H have consid¬ 
ered inflation driven by a tower of Kaluza-Klein fields with quadratic potentials. However, 
this mass distribution is very different from the Marcenko-Pastur form relevant in our case 
of string axions.) It turns out that while there are very general results for the perturba¬ 
tion spectrum produced by multi-field inflation ( e.g . P0fl), the spectrum of N-flation is 
surprisingly complex and has several previously unsuspected features. In this paper, we 
work only with the slow-roll approximation, and the results here represent an initial survey 
rather than a comprehensive analysis of the subject. We would like to emphasize that the 
dynamics and the resulting inflationary power spectrum do not become independent of the 
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initial conditions at late times [17|. Thus, even if one had a microscopically computable 
N-flation model, some of the predictions would depend on the initial conditions. For this 
reason, we present our results as functions of the initial conditions, illustrating the sim¬ 
plifications that arise in a few special cases. A theory of the initial conditions, which we 
do not address in this paper, would still be necessary for a complete understanding of this 
class of models. 

Readers whose main interests are the inflationary dynamics and the predictions of the 
model, rather than the problem of the axion mass spectrum in a string compactification, 
could begin with §2.1, examine the mass spectrum of §4.2, and then proceed from §6. 


2. The N-flation Model 

2.1. Basic Form of the Model 

Here we briefly recall the structure of the model proposed by Dimopoulos et ai |[I]| . 
Consider N axions whose periodic potentials 

y( W ) = A| (l-cos(^)) (2.1) 

are independent and arise solely from nonperturbative effects. As we will explain, we omit 
multi-instanton corrections to this potential. The scale A t is determined by dimensional 
transmutation and can be hierarchically small compared to the ultraviolet cutoff. We find 
it convenient to absorb the axion decay constants fi via 

4>i = -y, (2-2) 

Ji 

and we take units where M p — 1. For small axion displacements TC I , the potential 
may be Taylor expanded, and the Lagrangian isB 

L = (2.3) 

If the masses rrii are very nearly equal and the axions are initially displaced by a 
common amount (f>, they will roll in unison toward their joint origin. For all practical 
purposes, this is equivalent to the motion of a single field $ with a displacement y/Nf>. It 

6 For the moment we have neglected the important possibility of nontrivial axion kinetic terms, 
but we will come to a complete treatment of this point. 
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follows that even if each axion is displaced by a distance in field space small compared to 
the Planck mass, <f> can satisfy the conditions for slow-roll inflation. In this sense, N-flation 
realizes m 2 (j) 2 inflation in a very well-controlled string theory setting. 

One should ask how often this simple potential involving N decoupled axions with 
canonical kinetic terms and identical masses can arise in a string compactification. We 
would expect that a generic model with N decoupled axions will have a nontrivial spectrum 
of masses. Our goal is to determine exactly which sorts of mass spectra are possible in 
string compactifications. Surprisingly, we will find that a single characteristic spectrum 
emerges, regardless of microscopic details. 


2.2. N-flation in a Stabilized Vacuum 


An important requirement for a realistic inflationary model in string theory is the 
stabilization of compactification moduli. For this reason, it is important to verify in detail 
that the particular classes of axion potentials needed for N-flation do arise in a class of 
stabilized string compactifications. 

The moduli are typically paired with axions in complex combinations. In the type 
IIB theory, for example, the Kahler modulus p r measuring the volume of a given four-cycle 
B\ ; combines with the axion 


B 


(O 


Ca 


(2.4) 


as pi — i(j)i. (We have chosen a convenient sign.) 

There is a simple relationship between the moduli-stabilizing potentials and the axion 
potentials in the resulting stabilized vacuum. For definiteness, we will focus on moduli 
stabilization in type IIB string theory by the method of KKLT [Q, in which all moduli 
are stabilized by the combination of fluxes 0 and a nonperturbative superpotential. In 
the type IIA vacua of |31|], most fields receive masses from fluxes, but the lightest axions 
are lifted by nonperturbative effects; for this reason, we expect rather similar qualitative 
features in this context.0 We leave the properties of N-flation in more general stabilized 
vacua as an interesting subject for future work. 

The KKLT superpotential is 


W = W 0 (S,Xa) + YlMXa)e~ ai 

i 


7 This conclusion was reached in discussions with A. Krause and A. Mazumdar. 


(2.5) 
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Here 


Wo 


/ 


G 3 A O 


( 2 . 6 ) 


is the flux-induced superpotential, which depends on the dilaton S and on the h 2 ,i complex 
structure moduli Xa- The sum of nonperturbative superpotential terms runs over all the 
Kahler moduli Tj = pi — ifa (i = 1... hip). The real parts pi measure the volumes of 
four-cycles in the Calabi-Yau, and the imaginary parts d> % are the axions arising from the 
four-form potential C,\. We will work at volumes large enough to ensure that neglecting 
the multi-instanton corrections to ( |2.5| ) is consistent (see §9.1). 

We are assuming that there is one nonperturbative effect for each independent four¬ 
cycle. For each such cycle, there are two possible sources of a nonperturbative superpoten¬ 
tial: Euclidean D3-branes or strong infrared dynamics in a gauge theory, for which gaugino 
condensation on a stack of N > 1 D7-branes is the simplest example. The constants aq 
are either 2tc, in the case of Euclidean D3-branes, or 2tc/N for gaugino condensation. For 
simplicity we will assume that Euclidean D3-branes provide all the nonperturbative su¬ 
perpotential terms; relaxing this assumption is straightforward and does not affect our 
results. The prefactors Ai(xa), although knowable in principle, are unknown in practice. 
They stem from determinants of fluctuations on the D3-brane-instanton worldvolume. 

Finally, the Kahler potential is given, to leading order, by 



(2.7) 


where V is the Calabi-Yau volume, amounting to a highly nontrivial combination of the 


various Kahler moduli pi. (When there is only a single Kahler modulus, V cx p 3 / 2 .) 


The only feature of the KKLT scenario that will be important for our discussion is 
that the superpotential is a sum of an axion-independent term |32|] and N nonperturbative 
terms in which the axions appear as phases. This would still be true in an appropriate 
generalization of the KKLT method to other string theories. When this is the case, the 
axions potentials arise exclusively from the same nonperturbative effect that stabilizes the 
Kahler moduli with which they are paired. Our qualitative conclusions do not even depend 
on this basic property, but for the purpose of presenting a definite computation, we do 
assume the KKLT form except where noted. 
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2.3. F-flatness Conditions 

We are interested in a scenario in which the axions’ displacement from their joint 
minimum is the dominant source of positive energy during inflation. Once N-flation ends 
and the axions have relaxed to this minimum, supersymmetry may still be broken by some 
additional effect. However, this effect should be subleading, and ideally constant, during 
inflation, lest it ruin the useful properties of the axion-induced inflationary potential. 

The F-terms DaW of the various moduli S, y a , pi depend on the superpotential, and 
through it, on each of the axions. As the axions vary, the F-terms change, so the infla¬ 
tionary trajectory cannot be F-flat. However, we will substantially simplify our analysis 
by approximating the final vacuum as F-flat. This is an excellent approximation for phe¬ 
nomenologically viable models with low-energy supersymmetry: the F-term energy at late 
times can be at most of order the intermediate scale, and hence will be highly suppressed 
relative to the inflationary energy of the axionsJl Correspondingly, we will refer to config¬ 
urations as F-flat even if they satisfy only the weaker condition that the F-term energy is 
very small compared to the energy scales of interest for inflation. We could also imagine 
that the dominant supersymmetry-breaking in the final state is accomplished by a D-term 
effect such as an anti-D3-brane. In any event, we require that as the axions roll to their 
minima, the configuration relaxes to a KKLT vacuum, 

V(<l>f,S, X ,p) = V K KLT(S,x,p) ( 2 - 8 ) 

with (j)f the value the axions take at their minima. To summarize, then, our picture is that 
the axions are displaced from the minima they take in a supersymmetric AdS,\ vacuum, 
and their potentials are unaffected by any additional, subleading supersymmetry-breaking 
effects. 

We assume that all other moduli are fixed during inflation. By construction, all the 
moduli have nonzero masses in the KKLT vacuum, but we are requiring in addition that 
these masses are somewhat larger than the axion masses. A condition of this nature is 
essential to prevent unfixed moduli, whose potentials are not protected against corrections, 
from spoiling the properties of an inflationary epoch in which only the axions are dynamical. 
At first sight it may appear implausible that a modulus field pi would be firmly stabilized 
while its axion partner fa from the same chiral multiplet remains somewhat mobile. The 
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We thank S. Kachru for helpful discussions of this point. 



important difference is that the pi appear in the Kahler potential, whereas the (pi do not; as 
a result, when the KKLT potential is expressed in terms of canonically-normalized fields, it 
is steeper in the pi directions than in the (pi directions. Equivalently, the difference is that 
the (pi enjoy the protection of an approximate shift symmetry, but the pi , like the other 
moduli, do not. During inflation the unprotected moduli will typically acquire masses 
m> H, while the axion masses will remain substantially smaller.il 
To find the supersymmetric AdS 4 vacuum, we solve 

0 = D a W = d A W + (d A K) W (2.9) 

where A runs over the dilaton S, the complex structure moduli y„. and the Kahler moduli 
Pi. The F-terms ( |2.9| ) are axion-dependent, and cannot vanish unless <p = 0 or it. However, 
the F-terms also depend on the other moduli S, ;y, p. As explained above, while the axions 
are supporting inflation, the other moduli must remain fixed. Together with ( |2.8|) , this 
implies that the full axion potential is 

rw) = r(Apf, x f ) .s<")) (2.10) 

i.e. the supergravity potential evaluated on a background configuration in which pi,x a , S 
take the same values they take at the global minimum. Although the F-terms do not 
vanish along the axion trajectory, but only at the endpoints, one must still enforce 

D A W{(pi = 0) = 0 (2.11) 

This amounts to setting the axions to zero in each of the F-terms and using the residual 
relations to constrain Pi : Xa, S. These constraints will provide very useful clues about the 
structure of the axion potential. 


3. Mass Terms in N-flation 

We will now work out in some detail the potential arising from a KKLT realization of 
N-flation. 

9 See e.g. |k| for a discussion of related issues. 
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3.1. The Quadratic Potential 

Recall that the axion kinetic terms are identical in structure to those of the Kahler 
moduli with which they are paired, so that 

Caxion \M 2 v K l3 V^V'V - V (3.1) 

The supergravity potential is 

v = e M^)( KABD - WD * w - 3 ^f) < 3 - 2 ) 

where, as before, A. B run over the dilaton S. the complex structure moduli y„. and the 
Kahler moduli pi. (In this section we find it convenient to retain explicit factors of M p , 
which we have set to unity elsewhere.) 

In the KKLT superpotential ( |2.5| ), each nonperturbative factor depends on the asso¬ 
ciated axion as 

Wi = A ie - 27rpi e 2ni ^ = Cie 2ni ^ (3.3) 

We have chosen the notation Ci to stress the fact that the moduli y„, pi which appear in 
Ci are not dynamical variables in the scenario of interest: by ( |2.9| ) . these fields are fixed to 
the values they take in the supersymmetric minimum. Thus, the Ci are constants as far 
as the inflationary dynamics is concerned. Because the axion (p t appears only in W,. and 
neither in the Kahler potential nor elsewhere in the superpotential, every factor of e 2nl ^ i 
in ( |3.2|) will appear multiplied by the constant C r . 

Inserting ( |2 . 5|) in ( |3.2|) , performing a Taylor expansion around the origin (j)i = 0, and 
using the F-flatness conditions D^W |</, i= o = 0, we arrive at 

V = (2 ?TfMijptf + O(0 3 ) (3.4) 

where 

Mij - j^e K ( K AB D A CiD B C 3 - 3 C % C^ (3.5) 

The kinetic terms (El) and the mass matrix (|3.4|) are evidently not diagonal in gen¬ 
eral. Thus, we find that the potential does not automatically take on the uncoupled form 
required for N-flation. That is, the mass matrix is typically not diagonal in the the basis in 
which the superpotential is simple, with a single axion appearing in each nonperturbative 


10 






term.lij Moreover, the off-diagonal terms appear to be of the same order as the diagonal 
ones. We show in Appendix A that for essentially statistical (rather than directly physical) 
reasons, the cross-couplings in My are actually suppressed relative to the diagonal terms. 
Even so, keeping track of these cross-couplings will provide a key insight into the mass 
spectrum. 

At this stage, it is convenient to make a change of basis that renders the axion kinetic 
terms canonical. Define O to be the matrix that diagonalizes the metric K t ,j on the Kahler 
moduli space, 

OfKuO/ = k;X- (3.6) 


with k 2 the positive eigenvalues. The eigenvalues k, are related to the axion decay con- 
stants0 through 


Ki 


I 

M n 


(3.7) 


Next, rescale the <f) % to absorb the factors K t . so that the kinetic terms are canonically 
normalized and K tJ is the N x N identity matrix. Finally, redefine the complex structure 
moduli and the dilaton so that K ah , a,b — 1,... P is the P x P identity matrix, and return 


to units in which M p = 1. 


This leads to the Lagrangian density 


C 




(3.8) 


where 

M tj = (2n) 2 ^-O l k (D A C k D A C l - 3 C k C^C> 3 1 (3.9) 

Provided that the potential takes this purely quadratic form, i.e. neglecting the 
higher-order terms omitted in (|3]4|), there exists a further orthogonal rotation that diag¬ 
onalizes Mjj . The result is that all the data specifying an N-flation model can be stored 
in the list of eigenvalues of ( |3.9| ), which are the masses-squared rnf of N canonically- 
normalized, uncoupled axions. Our goal in the remainder of the paper is to characterize 
the statistical properties of these eigenvalues, based on the structure of (|3.9|). 


10 

li 


We thank S. Kachru and J. Wacker for useful discussions of this point. 


According to |34[], the Ki are at most of order unity; see also [p5|. 
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3.2. The Complete Potential 


In the previous section we focused on the quadratic potential (|I.4|). However, terms 
of higher order in the axions must be retained when some of the initial dimensionless 
axion displacements are of order one. When the quadratic approximation to the 
axion potential is invalid, it is generally impossible to diagonalize the potential to remove 
cross-coupling terms. 

Because large initial displacements are advantageous for inflation, there is a compe¬ 
tition between obtaining enough inflation and having a controllable expansion of the full 
potential (!T2) as a potential for uncoupled fields. For that reason, in this section we will 
present the full, unexpanded potential ( |3.2|) and sketch the resulting axion dynamics. 

To do this, we reorganize ( |3.2|) and use the F-flatness conditions ( |2.11| ). The result 
from the F-term piece@ of the potential is 


V — ^^(e K K ab DaCiDbC^] fl + cos(27r (pi — 2it(j)j) — 2 cos(27r (f)i) 


(3.10) 




Taylor expansion of this potential leads immediately to the F-term part of (|3.4|) . 
The potential ( |3.10| ) can be usefully rewritten in a more general form 

V = Vq + ^aiCos(27T0i) + ^ fJ t jCOs(2n(f>i - 2iuf)j) 
i i,j 


(3.11) 


from which we can recover ( |3.10| ) if we take the obvious definitions of the axion-independent 
term Vo and of the coefficients a L . 

The simplest sort of N-axion potential would have been a sum of cosine potentials, 
each involving a single field, with negligible cross-couplings, i.e. (3ij — 0. We have just 
seen that such a simple potential never governs deviations from an approximately F-flat 
KKLT minimum. While we do not claim that uncoupled N-axion potentials are impossible 
in string constructions, we do expect that they are rather rare in general. It would be very 
interesting to find exceptions to this rule. 

Requiring that the /3{j are negligibly small is also a very significant constraint on the 
general form ( |3.11 ). The problem is that there are N 2 cross-coupling terms but only N 
single-field terms, so cross-couplings which appear small could still be important in large 


12 The terms from —3|W | 2 have an analogous form but complicate the formulas substantially. 
Because these terms are unimportant for our considerations, we omit them here. 
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numbers. (However, as explained in [[1], the fact that each field appears in only a small 
fraction of the couplings suggests that assisted inflation is still possible. See the Appendices 
for a more detailed discussion of these issues.) 


In our particular case fl3.10|) , it is apparent that the net effect of the terms coupling two 
axions is not subleading to the effect of the terms involving a single axion. For example, 
when all the axions are near the origin (j)i — 0, the two classes contribute amounts of 
energy that differ only by a factor of two. If instead the axions have a common vev <f> ~ 
the cross-terms give the dominant contribution to the potential. Moreover, the energy 
from these cross-terms is positive and acts to repel the axions from the same point of their 
configuration space, substantially complicating the motion of the axions. The identification 
of a single direction for collective motion (such the inflaton p of [:!]) could be challenging. 


We will leave the full dynamics of the N-axion potential (|3.11|) as an interesting problem 
for the future. 

Finally, it is important to observe that diagonalizing the mass matrix is not equivalent 
to sending (3ij —> 0: the presence of cross-couplings ensures that the eigenvalues of M t] 
are different from the a t . The correct procedure, when the quadratic approximation is 
applicable, is to rotate away the cross-couplings, not to omit them. Surprisingly, this 
simplifies the problem of understanding the mass spectrum: we will see that although the 
collections {o. t } vary considerably from one compactification to another, the statistical 
properties of the spectrum of eigenvalues { rnf } of M %3 are much less variable. 


4. Random Matrix Treatment of the Mass Terms 

In order to understand which classes of N-flation models are most likely, we must 
understand the properties of the mass matrix M l3 , ( |3.9| ). However, the complete expression 
is prohibitively complicated, particularly because it depends on the unknown quantities 
DACi . For this reason, we resort to characterizing the expected properties of the ensemble 
Q m of mass matrices one expects in an ensemble of N-flation models. More specifically, we 
seek the probability distribution function p{m 2 ) for the eigenvalues mf of the mass matrix, 
a function which we call the mass spectrum for short. 

The logic behind this step is that, even without knowing the individual terms DACi , 
we can learn something about the eigenvalues of the mass matrix by searching for proper¬ 
ties of the distribution of eigenvalues of ( |3.iJ| ) that are independent of the particular set of 
values taken on by the DaCi in any given realization. This approach is effective because 
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the mass matrices are so large that fluctuations away from the mean behavior are unlikely. 
Thus, knowledge about average elements of Om provides an excellent guide to the prop¬ 
erties of nearly all the members of the ensemble. We will therefore restrict ourselves to 
characterizing Q m as a whole. 

Moreover, although we will be able to establish some facts about the statistics of the 
DaCi themselves, 0 it is possible, and advantageous, to find properties of the eigenvalue 
spectrum that are substantially independent not just of the particular values attained 
by the DaCi, but even of the distribution Od of these values in an ensemble of models. 
Because the matrices in question are moderately large, simplifications of the sort that often 
arise at large TV in random matrix theory turn out to simplify the problem. A particular 


result |27] which provides a strong motivation for our approach is this: the behavior of the 
eigenvalue spectra of TV x TV random matrices whose entries have appropriately bounded 
moments but an otherwise arbitrary distribution Qjj, is identical, in the TV —► oo limit, 
to the spectrum in the simplest case in which the entries are independent and identically 
distributed (henceforth i.i.d.) and is a Gaussian distribution with mean zero. This is 
an extremely powerful statement, and although we will not appeal to it directly, it gives 
an underlying reason that one might have expected universal behavior in our system. 


4-1- Statistical Model for the Covariant Derivatives 


In the basis in which the kinetic terms are canonical, the mass matrix takes the form 
( [1.01 ) . The mass terms arising from —3|!T| 2 are sub leading, by a factor of the inverse 
volume -, to those coming from the F-terms. For this reason, we neglect the former in 
the analytical discussions that follow. In §5.2 we will present the results of a Monte Carlo 
simulation of the full structure ( |3.9| ) to show that this approximation is a sound one. There 
we will find that even without including the suppression by A the terms from —3|IF| 2 have 
a negligible effect on the statistical properties of the mass matrix. 

The remaining quantities of interest are the rotation matrices OS, the decay constants 
fi, and, most importantly, H a% = D,\Wi- the (TV + P) x TV matrix of Kahler-covariant 
derivatives. 

We will now attempt to characterize the statistical properties of the entries in H. 
We propose that the entries of H may be modeled as independent, identically distributed 
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More sophisticated treatments of the statistics of flux vacua appear in [ 361. 
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(i.i.d.) variables with mean p// and variance a 2 u . We do not assume that the distribution 
is Gaussian. 

From these assumptions, one can show that in fact p# -C (?h■ Observe that the 
row-average ta = = — jjDaWq by the F-flatness condition. Granting the i.i.d. 

assumption stated above, the Central Limit Theorem implies that the row-averages r a are 
drawn from a distribution with mean /ih and standard deviation ^=. The sign of every 
ta is therefore the sign of fin , unless an > y/NH rr (so that a fluctuation can change the 
sign of one or more of the Ra •) 

However, genericity and microscopic considerations require that the terms r,\ — 
— jjDaWq take on both signs. In particular, the sign of = — jjWodiK depends on 
the intersection numbers CVp via d,K = ~ y QpVW't N-flation is only possible in com- 
pactifications in which some fraction of the C l3 k are negative (in order to avoid excessive 
renormalization of Newton’s constant H), so the ta cannot all have the same sign. If 
our distribution is to reproduce this basic feature, the mean value p# of the Hai must be 
taken to be parametrically small compared to the standard deviation, 


I^H 


< 


(?H 

Vn 


(4.1) 


For all practical purposes we may therefore set // // = 0.0 
Now define 

Rm = 2m K Vfr 1 0, i H Aj 


(4,2) 


so that the leading contribution to the mass matrix is 


Mij = R, iA R A j 


(4.3) 


Provided that the rotations O are random, it follows from //// = 0 that the mean // of the 
RAi may likewise be approximated by zero. However, the variance a 2 of the Rai may be 
different from a 2 u . 

In summary, we claim that for the purpose of determining the spectrum, the essential 
structure of the mass matrix is (|4.3|) , where Rai is the (IV + P) x N matrix defined 
in ( [4.2|) . We further propose to approximate the entries of R as i.i.d. variables with zero 
mean and some variance a 2 . We repeat that the entries do not necessarily have a Gaussian 
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Even if an exception to this result could be found, the addition of a mean yn ^ 0 would not 


change the spectrum of masses of the N — 1 lightest axions [27|. 
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distribution. To underscore both this point and the results of §4.1, we will construct, in §5, 


ensembles of random mass matrices using the full formula (|3.9|) . For a range of assumptions 
about the distributions of the entries of K AB and II a%- we will find results indistinguishable 
from those that follow from ( [4.3|) . 

We have not found an exact expression for a in terms of microscopic parameters, 
beyond the definition 

a={\RAi\) (4.4) 


where the average is over all values of the indices A, i. However, a useful heuristic is 

a2 ~W) {{DA Ci)2) (4 ' 5) 

which shows that a measures the typical size of the F-terms. The scale of the total super¬ 
potential, and thus of the F-terms, is determined by size of the flux-induced superpotential 
Wo, which in turn is fixed by a choice of fluxes. As in other applications of the KKLT 
proposal, a fine-tuning of the value of Wo is necessary in order to construct a vacuum at 
large volume and weak coupling. Although many moduli are present in our setup, a single 
fine-tuning, that of Wo, suffices to fix the overall scale, and hence the typical axion mass 
fh. In practical terms, we set the value fh exactly as in [[I], by requiring that the density 
perturbations resulting from inflation have amplitude consistent with observations. This 
implies (see §8) 

fh « 1.5 x 10 -5 M p (4.6) 


Let us point out that (|4.6| ) is not satisfied by an N-axion potential in a truly generic string 
compactification, which would presumably have an average mass of order one in string 
units. We must restrict our attention to the subclass of models (with fine-tuned fluxes) in 
which the overall scale does obey ( |4.6|) . Of course, this problem is not special to N-flation: 
in nearly every model of inflation in field theory or string theory, one must set the scale of 
inflation by hand. In our case, the necessary fine-tuning happens to be of the form treated 
by Bousso and Polchinski [^7 . 


The distribution of the axion masses around the average scale ( |4.6| ) is the quantity 
that we will now determine with random matrix techniques. 
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4-2. Axion Masses and the Marcenko-Pastur Law 


From the complicated structure of the mass matrix (|3.9|) we have abstracted one 
essential property: 

M = R t R (4.7) 

where R is a (N + P) x N rectangular matrix whose entries are i.i.d. with zero mean and 
variance a 2 . We have therefore reduced the problem of finding the mass spectrum in a 
typical N-flation model to that of finding the eigenvalue probability distribution function 
p(m 2 ) for a matrix of the form (|4.7|) . This important problem^ was solved by Marcenko 


and Pastur in 1967 [09 


The Marcenko-Pastur law for the eigenvalue spectrum is 

1 


p{m 2 ) = 


2 -7T m 2 /3a- 


-\J (b — m 2 )(m 2 — a) 


for a < m 2 < b. where we have defined 


(4.8) 


a = or 2 (l — l/tf) 

(4.9) 

b = CJ 2 (l + y/p) 

(4.10) 

N 

(4.11) 

y N + P 


The probability density vanishes outside this range. In Fig.l we display the Marcenko- 
Pastur law for a few illustrative cases. 


The parameters controlling the shape of the spectrum ( |4.8|) are A and a, where (3 
is the ratio of the dimensions of the rectangular matrix R and a 2 is the variance of the 
entries of Ra%- In the case of a N-flation in a KKLT compactification of type IIB string 
theory, there are hip axions, and hip + h 2 p + 1 is the total dimension of the moduli space 
(Kahler, complex structure, and dilaton), so that 


P 


h 


bi 


hi i + ho.i + 1 


(4.12) 


15 The problem of finding p( A) for an ensemble of matrices of the form (4.7) is important for 
statistical analysis and also has a range of practical applications, for example in communications 
engineering [ j38f . Here we denote the eigenvalue by A because the interpretation of an eigenvalue 
as a mass-squared is rare in these other applications. 
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Relative probability 



Fig. 1: The Marcenko-Pastur law for f3 — 9/10 (highest peak), 1/3, 1/10 
(most localized). 

In more general cases, for example in other string theories, it may be possible to identify 
f3 with a ratio of this form: the number of (real) axions involved in inflation divided by 
the total number of moduli chiral multiplets. 

We will not attempt to find the smallest and largest values of (3 that can arise in 
a working N-flation model. In any case, there is an important, independent constraint- 
on (3. The renormalization of Newton’s constant described in [;1|] and reviewed in §9.1 
requires that N ~ P, so that (3 ~ 4. In special (and perhaps fine-tuned) circumstances 
this constraint might be relaxed, so we perform the remainder of our analysis for a general 
(3. However, we will keep in mind that models with (3 ~ \ are strongly favored because 
they more readily produce an appropriate number of e-folds. 

To determine the value of cr, notice from ( |4.8|) that the average value of the mass- 
squared is 

(m 2 ) — a 2 (4.13) 

from which it is evident that a sets the overall mass scale. As explained in the previ¬ 
ous section, the only phenomenologically viable N-flation models are those in which ( f4.6| ) 
is satisfied. This effectively fixes the value of a and leaves f3 as the unique parameter 
specifying an N-flation model. 
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4-3. Some Relations to Other Results in Random Matrix Theory 


The Marcenko-Pastur law [39] is in some sense analogous to Wigner’s celebrated 


semicircle law |25[ for the distribution of eigenvalues of a random symmetric matrix 
M — H + H T , where the entries of the square matrix H are i.i.d. Gaussian variables 
with mean zero. In the same way, Q m has some superficial similarities to the famous 
Gaussian Orthogonal Ensemble (GOE), which describes symmetric matrices whose entries 
are drawn from normal distributions. 

Recall that elements of the GOE []24|] may be constructed as G — H + II 7 , where 
the entries of H are i.i.d. Gaussian variables with zero mean. The eigenvalues of G obey 
Wigner’s semicircle law. Our mass matrix has instead the structure M — II 1 II. and 
substantial correlations among the matrix entries in M conspire to give the Marcenko- 
Pastur spectrum of eigenvalues. 

The Marcenko-Pastur law (|4.8D does approach a semicircle law, although not one 
centered at the origin, in the limit (3 —> 0. This limit corresponds to P/N —> oo, so that the 
number of axions is very small compared to the total number of moduli. Compactification 
manifolds in which P 3> N and N 1 are presumably quite rare, so we find that the 
Marcenko-Pastur law is quite different from the semicircle law for realistic values of (3. 

In the other limiting case (3 —> 1, the Marcenko-Pastur spectrum approaches the 
square of a Wigner semicircle spectrum (also known as the quarter-circle law |j2G| ); that 
is, one can construct the former spectrum by squaring the eigenvalues appearing in the 
latter. Squaring a symmetric but otherwise general matrix H is one way to construct a 
matrix which is positive-definite (as we have demanded that our mass matrix should be) 
and has an easily-computed spectrum. In the limit (3 —> 1, in which there are very few 
chiral multiplets in the theory except for those in which the axions appear, the N-flation 
mass spectrum actually arises in this very simple way. 

In the special case that the distribution of the Ra% is Gaussian, M is actually a real 
Wishart matrix [ID|, about which a great deal is known. In this paper we do not assume 
that the distribution is Gaussian, as there is no clear physical motivation for this beyond 
the universal behavior of random matrices 10 in the N —> oo limit. However, in the event 
that exact results are desired for a toy model, the Wishart ensemble could prove to be 
quite tractable. 
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5. Monte Carlo Analysis of Random Mass Matrices 

In the preceding section we made a few assumptions about the statistical properties 
of the Kahler-covariant, derivatives D,\C r and about the Kahler metric itself. Moreover, 
we argued that the essential structure of the mass matrix (|3J^) is captured by (|4.7|). These 
simplifications led to a universal result for the mass spectrum, the Marcenko-Pastur law 


In this section we use Monte Carlo methods to provide further evidence for the con¬ 
clusions of §4. In particular, we show that the results are robust against relaxation of our 
assumptions. We construct a collection of mass matrices according to the full rule (|3.5|), 
using a variety of input distributions for DaCi,K ab . The result is indistinguishable from 
( |4.8|) , providing a strong confirmation our analytic result. 

To generate a random mass matrix M r ,j . we choose the entries of the component ma¬ 
trices K ab , CiCj , and D,\C, from appropriate ensembles Ak- 0^', - hi- This procedure is 
clearly superior to postulating the distribution Om of the entries of A/ iy directly, because 
the latter procedure produces no new information. In contrast, assembling the component 
matrices according to ( |3.5| ) may be expected to yield results that are substantially inde¬ 
pendent of the (largely unknown) details of Aq- -hr- but that do depend on the known 
structure of the matrix products in ( |3.5|) . We will see that this expectation is justified. 


5.1. Properties of the Component Matrices 

Although the Kahler metric is determined by the intersection numbers CVp : of the 
Calabi-Yau, the scarcity of examples where all intersection numbers are known makes 
this an impractical route to an understanding of the properties of a typical K IJ . For this 
reason, our examples of “Kahler metrics” were simply random, symmetric, positive-definite 
matrices. 

The derivatives D^Ci split into two classes. For A — a, i.e. when the derivative is 
with respect to a complex structure modulus, we have 

D a C i = d a KC i + d a C i (5.1) 


and we recall that Ci may depend on the complex structure sector through the prefactor 
A, which is a fluctuation determinant in the case that the nonperturbative superpotential 
comes from Euclidean D3-branes. As we have stressed, the form of this dependence for a 
general threefold is not known. 
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For A — i, i.e. in the Kahler sector, we have instead 


D,,Cj — diKCj — 2nCj8ij (5.2) 

so that the average behavior of the diagonal elements is distinguishable from that of the 
off-diagonal elements. 

Observing that each of the three classes of individual terms d a Ci, BaK, Ci can be 
positive or negative, we constructed samples of DaCi by drawing these terms from dis¬ 
tributions that ranged over both positive and negative numbers of order one. We found 
that for a wide range of properties of these distributions, the mean value /i of the elements 
of the resulting matrices DaCi is much smaller than the standard deviation a. This is 
consistent with our analytic demonstration, in §4.1, that the mean value of these elements 
is necessarily small compared to their standard deviation. 


5.2. Results of Monte Carlo Simulation 


We generated thousands of examples of mass matrices using the rules outlined in the 
previous section. In particular, we varied the distributions Qk,^c,^h- For example, in 
the case of Qh we considered normal distributions with various values of the mean; uniform 
distributions; log-normal distributions; chi and chi-square distributions; and many other 
distributions assembled from these via matrix products. 

The results of this simulation are as follows. The distribution of matrix elements is 
bimodal, with the off-diagonal entries normally distributed around zero, and the diagonal 
entries normally distributed around a positive mean. As may be expected from the central 
limit theorem, the simulated distributions are extremely close to normal distributions for 
a wide range of input parameters. 

The positive mean of the diagonal terms is large compared to the expected magnitude 
of the off-diagonal terms. Specifically, we find that the typical magnitude of off-diagonal 
elements is smaller than the typical magnitude of diagonal elements by a factor 


1 

VnTp 


(5.3) 


The numerical coefficient is of order one. This suppression is of limited relevance for our 
primary goal of understanding the mass spectrum, but it could be useful in other contexts, 
such as an analysis of the full axion potential ( |3.10| ) . For this reason, we provide an analytic 
derivation of this suppression factor in Appendix A. 
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Finally, and most importantly, the eigenvalues are indeed distributed according to 
the Marcenko-Pastur law, as shown in Fig. 2. The shapes of the empirical eigenvalue 
distributions agree extremely well with those predicted by (|4.8|) , without any need for 
tuning or fitting of parameters. 

The spectra shown are for the simplest case, that of a Gaussian distribution with mean 
zero. In each of the other cases we tested, the spectra are identical to those shown, except 
for the presence of a single eigenvalue N times larger than the others. The presence of 
this single large eigenvalue is familiar from random matrix theory [^7|. In physical terms, 
the associated field will rapidly move to its minimum and decouple from the inflationary 
dynamics. Because a spectrum with an extreme outlier is not easily interpreted visually, 
we do not show an example. 


6. N-flation Dynamics and Initial Conditions 

Having derived the mass spectrum for the fields contributing to N-flation, we are now 
in a position to evaluate the cosmological evolution associated with this system. 


6.1. Inflationary Dynamics 

As usual, the dynamics are described by the general set of equations, 


H 2 = 


1 


3 Ml 


+ v ((p) 


2 


dV 

+ 3H(pi + —— = 0 

ofa 


( 6 . 1 ) 


( 6 . 2 ) 


(6.3) 


We will make frequent use of the fact that our potential V is written as a sum over terms 
of the form m 2 (f) 2 /2 . 

During slow roll, one can find a general result for the held evolution [|17|. Dropping 
the acceleration terms from the fields’ equations of motion, and dividing the fa equation 
by the (f)j equation, we see 

fa rrif fa 


°J "> j °i 


(6.4) 
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Fig. 2: Empirical eigenvalue spectra of 100 random 300 x 300 matrices 
constructed according to ( |4.7| ), with (3 — 1/10, 1/3, 9/10. The horizontal 
axis is the mass-squared and the vertical axis is the relative probability. The 
overlaid curves show the Marcenko-Pastur law in each case. 
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*?A 


Mto) 

This is not an attractor: if we rescale the initial value of one of the fields by multiplying 
(f>i(to) by a constant, this rescaling changes the held value at all subsequent times, for as 
long as the slow roll approximation remains valid. 

It is instructive to specialize to the Marcenko-Pastur spectrum, where the lightest 
held 4> 1 has m 2 ~ a (ED, the average mass-squared is a 2 , and the most massive held (J)n 


has m 2 N ~ b (|4. 10|) . Let us use the progress of (J)n as a clock, via 


T(t) = 


4>N(to ) 


( 6 . 6 ) 


As inflation proceeds, r decreases. Trivially rewriting (|6.5|) , we have 




(6.7) 


where Q % = -r- are in the range 


a 

b 


1 + \//7 


< Q < 1 


( 6 . 8 ) 


(In the preferred case (5 ~ we have f ~ -yj-) Even if the helds begin near the same 
point in held space, they quickly spread out: the large range in the exponents Qi ensures 
that the helds move at very different rates. We illustrate this spreading in Fig. 3. 

If instead each held had a potential, the solution analogous to ( |6.5|) would be 


1 _ 1 _ Xi_ f 1 _ 1 

( kit ) 2 (j)i{t 0 ) 2 A j \(pj(t ) 2 (pj(t 0 ) 2 

We observe that the initial conditions become increasingly irrelevant at late times, demon¬ 
strating the existence of an attractor in the Ap 4 case. This is a desirable feature in a 
model of inflation, but because realizing assisted Ad 4 inflation in string theory (or in any 
well-motivated particle physics setup) is probably very difficult, we will say no more about 
this case. 

The absence of an attractor solution for N-hation signihcantly complicates the inha¬ 
tionary dynamics. If there were a late-time attractor, and inflation lasted long enough to 
eliminate the residual effects of the initial conditions, we could assume that the system 
followed the attractor solution during the last 60 e-folds of inflation, greatly simplifying 
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Fig. 3 : We plot the evolution of fa as a function of time for N — 201, 
j3 — 1/2, and i — 1,21, •••,201 with m — 10 _6 Mp, and time measured in 
units where m — 1. 

the treatment of the perturbations. By contrast, we will see in the next section that the 
N-flation spectrum is a function of the values of the fields measured at some fiducial time, 
since the initial conditions never become irrelevant. 

Finally, notice that in single-field inflation, the breakdown of slow roll coincides with 
the end of inflation, and to a first approximation the two events can be regarded as equiv¬ 
alent. In N-flation, however, slow roll can fail well before the end of inflation for fields that 
are particularly massive or that begin close to the origin. 

6.2. Initial Conditions 

We will not give a full treatment of the initialEl conditions in this paper. However, a 
few observations will be essential. 

The data of an N-flation model in the four-dimensional effective theory is a spectrum 
of N axion masses. Given any such system, one would like to know whether the model can 
produce, with some choice of initial conditions, J\f > 60 e-folds of inflation. For this purpose 
it is useful to choose initial conditions that maximize the number of e-folds while remaining 
consistent with all control assumptions, especially the requirement that the largest vev is 

16 Throughout this paper, we use “initial” to refer to the situation 60 e-folds before the end of 
inflation. This is not necessarily the start of cosmic history. 
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at most of order M p . The amount of inflation is maximized by placing all of the fields 
as far from the origin as this condition will permit; spreading the initial vevs inevitably 
moves some fields closer to the origin and thus diminishes the number of e-folds. N-flation 
systems do not readily produce a vast number of e-folds, so only a tiny fraction of N-flation 
models can generate enough inflation from a highly “inefficient” initial configuration with 
disparate vevs. For this reason, we will focus primarily on initial conditions in which the 
fields are spread by at most an order of magnitude. 

Consider any distribution {4>i} of initial field values, and compute the correlation 


_ (m 2 4> 2 ) _ NY, l m 2 4> 2 

~ (m 2 )(ct> 2 ) ~ Ej Efc 


( 6 . 10 ) 


where we have defined 

(m p 0 9 ) = -^r (6.11) 

i 

When £ > 1 (£ < 1), we say that the vevs are directly (inversely) correlated with the 
masses. 

We may also parametrize the correlation between masses and vevs as 


4>i = Vim? 


( 6 . 12 ) 


where the variation in prefactors v is presumed to be less important than the dependence 
on the mass. Let us try to relate a to £. We begin by observing that for reasonable values 
of N, (rn 2k ) is very well approximated as a moment of the Marcenko-Pastur distribution, 


(m 2k ) ~ TO( 2 fc) = / m 2k p(m 2 ) 


2 fc„ 


(6.13) 


where the approximation omits terms of order -jr 0 , which we can safely neglect in the 
following. The moments of ( |4.8| ) are in turn given by a simple analytic formula ||42|| , 


OU 

m {2 k) = 



(6.14) 


We will also need to know (m 2 ) and (m 4 ). These are not moments of the distribution, 
but we can compute them directly from the Marcenko-Pastur law, finding 


(m 2 ) 


1 1 

^0~T) 


(6.15) 
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and 



(6.18) 


Equipped with these results, we can compute £(ct) in terms of (3 for the cases of interest. 

As a consequence of ( |6.5| ). the fields rapidly spread apart (except when the masses are 
degenerate), with the most massive fields relaxing most quickly. This means that the lighter 
fields soon have larger vevs, on average, than the heavier fields. Thus, positive correlations 
£ > 1 are transient, and inverse correlations £ < 1 typically develop. This situation is quite 
plausible dynamically, but, within the context of controllable N-flation models, it creates 
a certain tension. Initial conditions with £ > 1 are the most efficient at producing many 
e-folds within the region of control, 0 but because such positive correlations are typically 
transient, they may not be the most generic initial conditions. We will work not with the 
initial conditions that are most plausible a priori , but with those that are most plausible 
given the requirement of J\f > 60 e-folds. 

These considerations lead us to focus on initial conditions with a ~ 0 £ ~ 1, 

because anything else will typically give inadequate inflation. We will frequently refer 
to this as an “uncorrelated” initial condition, because the vevs are not correlated to the 
masses. The equal-field condition f>i — £> is an important special case of uncorrelated 
initial conditions. Another instructive configuration, albeit an impractical one for realistic 
models, is a — —1, where each axion makes a roughly equal contribution to the initial 
energy. 

We have numerically evaluated the full equations of motion for a wide variety of 
parameter values and initial conditions. In Fig. 4 and Fig. 5 we show the results of the 
evolution for two specific choices of initial conditions: the “equal-field” condition a — 0, 
v % ~ v. and the “equal-energy” condition a — —1, ~ v. 

6.3. Number of E-folds 

We will now study how the amount of inflation depends on f3. Armed with the general 
result & for the field evolution, we can write the potential V(t) as an appropriately 
weighted integral over the Marcenko-Pastur mass distribution, since 



(6.17) 


17 It might be that inflation occurs in nature in a way that is not theoretically computable; for 
obvious reasons we do not consider this case. 
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Field Evolution 
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Energy of Field / Dominant Field 



150 


Actual Evolution/Slow Roll Solution 



150 


N 



Fig. 4: N-flation dynamics with N — 200, (3 = 1, and m — 10 _6 Mp, 
and equal-field initial conditions. Time is measured in units where fh = 1. 
The fields are arranged along the axis labeled i, in order of increasing mass. 
The top left plot shows the field evolution. The top right plot displays the 
relative contribution of each field to the overall energy density, normalized 
by the field that is currently making the dominant contribution. The bottom 
left plot shows the ratio of the numerical solution to the slow roll solution, 
( |6.5| ), computed in terms of the vev of the lightest field. The bottom right 
plot shows the evolution of the scale factor as a function of time. 
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Field Evolution 



200 


Energy of Field / Dominant Field 



200 



Fig. 5: N-flation dynamics with N — 200, (3 = i, and m — 10 _6 Mp, 
and equal-energy (ct = —1) initial conditions. The plots follow the layout of 
Figure 4. This is an illustration, not a controllable model: in order for the 
total number of e-folds to match Figure 4, the initial vevs of the light fields 
have been taken to exceed Mp. This is a concrete example of the general 
observation that for a given mass distribution, initial conditions with a < 0 
require larger maximum vevs in order to produce enough inflation. 
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This is approximated by 

v ( r ( t )) = Y J Ax A p(A) (6.18) 

where 

x = r{t) 2 t b (6.19) 

and we have assumed that (j) t (to) = 0 for all i. However, this integral cannot be done in 
closed form, so we need to make a further approximation in order to compute the number of 
e-folds of inflation we can expect for a given initial value <p. In practice, we are particularly 
interested in moderate values of (3 - as (5 approaches unity, the mass distribution develops 
a very long tail, but this situation is disfavored physically. Consequently, we can get a 
rough estimate of the number of e-folds by assuming that the N fields have equal masses 
fh 2 — a 2 and possibly different initial vevs (pi . In this case, one can define <F 2 = N (<p 2 ) 
and write down an equivalent single-field model [jl] with potential 

H($) = -fh 2 $ 2 (6.20) 

2 


The power of assisted inflation is that each of the <pi is related to $ by <pi « $/ \/~N. and 
while the collective field $ must exceed Mp by at least an order of magnitude to obtain 
a workable period of inflation, the vevs of the <pi can all be sub-Planckian if N is large 
enough. 

The total number of e-folds in this simple model is 


N{$) 


i r $ 2 

mJ [t 


M p] _ $2 

IT] ~ 4Mp 


( 6 . 21 ) 


This turns out to be an excellent approximation to the number of e-folds produced in 
the more general case of a spectrum of masses. The quality of this approximation is 
demonstrated in Figure 6; we conclude that it is safe to ignore /3 when estimating A/", unless 
/3 is very close to unity, and provided that we take initial conditions with approximately 
equal fields. 

The number of e-folds also depends on a. As we are physically constrained to consider 
only sub-Planckian vevs for the (pi we can compute the maximum number of e-folds by 
saturating this bound for all the fields. Then any physically realistic initial configuration 
will have N vevs which are less than or equal to their values in this idealized scenario, and 
the amount of inflation produced when some vevs are initially reduced from their maximal 


30 







e-folds 



Fig. 6: Numerical results for the number of e-folds J\f as a function of f3, 
with N — 300 and initial field values cy — Mp. The expected number of 


e-folds from the approximation in ( |6.21| ) is 75, so when p is not very close to 
unity, the approximation is excellent. 


values will necessarily be lower. We can therefore find an upper bound on J\f for a given 
mass distribution by assuming that a- = 0 and (p l ~ Mp. 

Note that the mass scale fh = a has fallen out of ( 6 . 21 ). In the end, the mass is fixed 
not by the number of e-folds, but by the amplitude of the perturbation spectrum, which 
we now discuss. 


7. Generation of Adiabatic Perturbations 


There is a very general result, derived in 
bations generated by multi-field inflation, 


30|1 , for the spectrum of curvature pertur- 


Pn 


(H\ 2 8UdU 

\ 2it / dpi dpj l> 


(7.1) 


where J\f counts the number of e-folds. We will follow the treatment in 
izes ( |7.1|) to the case of uncoupled fields with polynomial potentials, 


I0|1, which special- 


r «■) = !>(*) 

i 


(7.2) 


In the case of N-flation, Vi — 
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The spectrum is given by 


101 


pi/2 

1 iz 


V 


N 


127r 2 Mf, 



(7.3) 


We have chosen to work with P,^ 2 , s i nce this is proportional to the amplitude of the 
primordial density perturbations, and to the amplitude of the temperature perturbations in 
the microwave background. We are implicitly ignoring entropy (isocurvature) fluctuations 
in this initial analysis, although these are likely to be subdominant for uncoupled fields 
with very similar decay channels. 

For our specific potential, the spectrum becomes 


pl/2 _ 
1 K — 



Combining these results with ( |6.11|) , we find 


(7.4) 


p i /2 = I (m 2 (j) 2 ) (cj) 2 ) 

n Mp V 96t r 2 


(7.5) 


We stress that ( |6.11|) is only a convenient notation, and does not imply that an aver¬ 
aging has taken place. The expression ( |7.5| ) is exactly equal to ( |7.4| ), for any collection of 
masses and field vevs. Of course, our purpose in writing (|775|) is to make use of statisti¬ 
cal information about the masses and initial conditions. To accomplish this, we will now 


use the moments ( 3.14 ) of the Marcenko-Pastur distribution, in combination with various 
possibilities for the initial conditions. 

In the case £ = 1 -<=>- a — 0 that the initial field values are uncorrelated with the 
masses, which includes the important special case of equal initial vevs, we have 


(m p (f) q ) = ( m p )((f) q ) 


(7.6) 


so that 


1/2 _ N^ 2 )^ 2 ) 


pi/* _ 


7T 


\/96Mj 


(7.7) 


If the masses are taken to be identical, this reduces to the result of an effective single-field 
model, 

,1/2 _ $ 2 to<i> _ V($) 


P ' = 

7rv^96M^ m$7T\/24M^ 


(7.8) 


32 


























pl/2 



• 1 / 2 

Fig. 7: Scalar power spectra P^ , computed from ( |7.4f) for the initial 

conditions used in Fig 4 (lower line) and in Fig 5. The masses are identical 
in the two cases. 


once we identify <F 2 = N((j) 2 ). This is the simplest possible model of N-flation. 

Considering instead the case a = — 1 in which each field makes a roughly equal 
contribution to the energy, we find instead 


pl/2 

1 n 


_ V _ 1 

mny/VAM^ y/1 — f3 


(7.9) 


which is enhanced relative to the single-field result by the modest factor (1 — (5) 1 1 2 r^j y/2. 

To get a clearer understanding of the connection between the scale-dependence in the 
spectrum and the initial field configuration, we now turn to the spectral index n s , which 
is defined as n s — 1 = dlnP-jz/dlnk, so that Pn(k) ~ k n3 ~ l . For a single inflaton, one has 
as usual 

n s — 1 — 2r/ — 6e (7.10) 


where 



In the case of single-field m 2 q t inflation, this becomes 


n, - 1 = 


8 Mi 


(7.11) 


(7.12) 


For multiple fields, n s obeys the general formula |[Tj 


n. 


1 = 2 


H 

JP 


\r ( i t\/t 2 Vij \ , r 




(7.13) 
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Using once again the specialized results of [10]lin this reduces to 


n s - 1 = 





(7.14) 


which reverts to ( |7.12|) in the single-field limit. For the special case of multiple quadratic 


potentials, the final term is proportional to the second term, so 


n s — 1 




(7.15) 


It follows that n s — 1 is never positive, and thus that the spectrum is never blue, because 
we have written n s — 1 as a combination of squares with a negative overall coefficient. 


n s 



e-folds 


Fig. 8: The spectral indices n s corresponding to the spectra plotted in Fig 
7, displayed as functions of log a. The case with equal initial vevs initially 
has a higher index. 


Using (|6.11|) we can rewrite (|7.15|) 


as 


n s — 1 = 


4 Ml 


1 (m 4 4> 2 ) 


N \(4> 2 ) ( m 2 (j) 2 y 


(7.16) 


Again, we stress that ( [7.16|) follows identically from ( |7.15|) , rather than being a statisi- 
cal approximation. Results that do use statistical reasoning are easily noticed because 


18 


This simplification is only possible when all the fields are still rolling slowly. 
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they contain factors of /?, the parameter that encapsulates the form of a typical mass 
distribution. 


For the case a — 0 of initial conditions that are uncorrelated with the field masses, 


we use (|7.6|) to find 


n. 


1 = 


8 Mi 




(7.17) 


N(4> 2 ) 

This coincides with ( |7.12| ) only in the degenerate limit j3 —> 0. Thus, for this important 
class of initial conditions, the generalized N-flation spectrum is more red (i.e. further from 
scale-invariant) than the single-field spectrum. This is one of our main results. 

In Fig. 9 we show the /^-dependence of the spectrum ( |7.17| ) in the case a — 0. Because 
the result shown is the spectral index more than 60 e-folds before the end of inflation, 
the values shown are closer to scale-invariant, than in observationally-relevant, cases. Our 
purpose in showing the figure is establishing the perfect agreement between ( 7.17 ) and the 
result of the numerical evaluation of ( |7.15|) . 


n s 



Fig. 9: The initial value of the spectral index as a function of /?, with all 
N — 300 fields initially equal to Mp. The solid line is the theoretical curve 
( FTP , and the dots come from numerical evaluations of ( |7.15|) . 


The quantitative relationship between n s and (3 depends on the initial conditions, and 
in particular on the correlation a between the masses and the initial vevs. As we stressed 
in §6.2, the displacements might become correlated inversely with the masses, because the 
more massive fields relax more rapidly. It would be very interesting to undertake a careful 
study of the initial conditions for N-flation, and to extract n s as a function of a and f3. 
However, this is a problem beyond the scope of the present paper. 
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Finally, we can compute the running dn s /d\nk of the scalar power spectrum, by 
applying the above methods to the general formula given in m ■ The leading contribution 
is 


dn s 

d In k 


16 

'W 


(mrcfr) 2 ( nrrcp 2 ) 


+ 


1 


(m 2 (f 2 ) 4 (m 2 </> 2 ) 3 (cf) 2 ) 


(7.18) 


For the uncorrelated initial condition a — 0, this simplifies to 

dll* 


dink N 2 (4 > 2 ) 


16M 4 9 

p (2 + f3 + f3 2 ) 


(7.19) 


This reproduces the single-field result if we again identify <f> 2 = N(cj) 2 ) and send f3 —> 0. 
Evidently the Marcenko-Pastur spread in masses creates a modest increase in the running, 
by a factor that is at most 2 and is typically around 

This analysis is far from being the last word in understanding the inflationary phe¬ 
nomenology of N-flation. To fully understand the spectrum arising from an arbitrary set 
of initial conditions, it will be necessary to go beyond first order in slow roll - and perhaps 
to numerically solve the evolution equations for the field perturbations, since the heavy 
fields are not well described by the slow-roll approximation during the last e-folds of in¬ 
flation. Moreover, because there is no attractor solution and the spectrum is a function 
of the initial field configuration, a solid prediction for the spectral index will require an 
understanding of the expected distribution of initial values for the axion fields. A further 
complication we do not consider is that the value of H at the end of inflation depends both 
on f3 and on the initial values of the axion vevs. In general, increasing (3 while holding 
the vevs and the overall mass scale fixed will lower the value of H at the end of inflation. 
This effect must be incorporated in a full calculation of the perturbation spectrum, as it 
changes the number of e-folds before the end of inflation at which a given comoving scale 
leaves the horizon, modifying the predicted spectral amplitude for that mode. Finally, it 
will be important to understand the possible role of isocurvature fluctuations. 


8. Astrophysical Constraints on the N-flation Parameters 

The simplest realization of N-flation is that in which N uncoupled, canonical fields 
have identical masses. As explained in [ill], this model is equivalent to chaotic inflation with 
the potential m 2 (j) 2 . On the other hand, we have also seen that a spread in the masses will 
modify the scale-dependence of the spectrum. We now turn to the observational constraints 
that can be imposed on the parameters of N-flation. Since the spectrum depends on the 
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initial field distribution, we will focus on qualitative constraints, rather than constructing 
detailed exclusion diagrams. 

First, we can normalize the mass scale by comparing Pp to the amplitude of the 
temperature fluctuations observed in the CMB 


Pr = 


800tt 2 
1 CMB 


A(k 0 ) = 2.95 x 10" 9 xl(fco) 


( 8 , 1 ) 


where A(ko) parametrizes the height of the spectrum at some fiducial wavenumber A: 0 and 
Tcmb is expressed in micro-Kelvins. Observationally, WMAP measures A « 0.78 H.F5 


For any given initial conditions the mass m = a is determined by comparing (7J3) and 
(0); however, the precise result can change by a factor of a few, depending on a and j3. 
Taking A/" ~ 60 and initial conditions that are not correlated with the masses (a — 0), we 
find 

1.5 x 10“ 5 M p (8.2) 


a 


In addition to the scalar perturbations, we must also consider the production of tensor 
perturbations during inflation. These are generated by quantum-mechanical fluctuations 
in the spacetime background itself, and have amplitude proportional to H cx W. One 
can seek these at CMB scales, through their contribution to the temperature anisotropy, 
through their contribution to the polarization B-mode (see e.g. (46|] for the likely sensitivity 
of future experiments), or via a direct detection of a stochastic background of gravitational 
waves with a mission such as BBO (47]]. Like the single-field m 2 (j ) 2 model it mimics, In¬ 
flation occurs at a comparatively high scale, and thus yields a significant tensor spectrum. 
However, to break the degeneracy between N-flation and single-field rn 2 © 2 inflation, a 
very accurate measurement of the tensor spectrum would be required. The exact relation 
between the tensor perturbations in the single-field model and in N-flation will depend on 
the distribution of initial vevs, as well as on j3. In practice, however, the amplitudes are 
certainly of the same order in the two cases, though not necessarily identical. Moreover, as 
both N-flation and m 2 (j) 2 produce a large tensor spectrum relative to most other inflationary 
models, a sufficiently tight constraint on the tensor amplitude would rule out both N-flation 
and single-field m 2 (j ) 2 models. 

The next constraint we explore is the lower bound on N, arising from ( |6.21|) and 
from the requirement that the individual axion vevs fa never approach Mp. Again taking 
M ~ 60, we find 

© 2 9dn 

(8.3) 


M 2 
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at the onset of inflation. Unless one can find reasonable string theory constructions that 
allow values of N as high as 10, the axion vevs in a model that produces adequate inflation 
will have to begin uncomfortably close to the Planck scale, though still below it. Although 
this situation is not entirely satisfying, it is still preferable to a single-held model in which 
the initial vev is parametrically large compared to the Planck scale, as it is improbable 
that models of the latter sort can exist in any controlled framework in string theory. 

In the previous section, we showed that the spectral index is similar to that produced 
by single-held m 2 (j) 2 inflation, but with somewhat more tilt at CMB scales (at least for 
those initial held conhgurations we considered explicitly.) These changes in the scalar 
power spectrum can be noticeable in the theoretically-favored case of the Marcenko-Pastur 
distribution with (3 ~ in that n s — 1 can change by around 25%. Observing this difference 
unambiguously would require a CMB / Large Scale Structure dataset of exquisite precision, 
but is potentially within the reach of anticipated CMB missions. 

In addition to curvature perturbations, multi-held inhationary models can also gener¬ 
ate isocurvature (entropy) perturbations,El as well as non-Gaussianities |49|| beyond those 
found in the single-held case. In this initial analysis we have ignored both of these effects, 
along with the possible impact of higher-order corrections to the slow roll formalism. Tak¬ 
ing these issues in reverse order, working at lowest order in slow roll gives us a very good 
account of the inhationary phenomenology of the single-held m 2 (j) 2 model. In N-hation the 
spread in the mass spectrum along with a generic, non-identical, set of initial vevs for the 
axions means that some helds will cease to be critically damped before the inhationary 
era comes to an end, i.e. before the second derivative of the scale factor a(t ) becomes 
non-positive. If the initial vevs are similar and (3 is not very close to unity, this breakdown 
in slow roll will likely only occur near the end of inflation, but in the general case one 
needs to exercise caution when using any slow-roll results. Further, understanding the 
non-Gaussianities [ |40| and isocurvature modes |5(J generated during multi-held inflation 
requires a calculation that extends well beyond the slow-roll approximation of the curva¬ 
ture perturbation, but may provide further avenues for breaking the degeneracy between 
N-hation and its single-held limit. We intend to address this problem directly in future 
work. In practice, a definitive understanding of the perturbation spectrum might require 


19 See [ ^8| for an interesting treatment of isocurvature in an assisted chaotic inflation model 
with O(N) symmetry. This analysis corresponds to the /3 —* 0 limit in our terminology. With 
(3 1/2 the masses vary significantly and the O(N) symmetry is badly broken. 
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evolving the N fields’ perturbation equations alongside the background equations, as was 
done in [[TT] for the single-field case. This would not be prohibitively challenging: the 
iV-field simulations displayed in this paper take a few seconds to complete, and adding the 
perturbations would roughly double the number of relevant degrees of freedom. Thus, the 
generic numerical problem may be computationally tractable. 

Finally, it is encouraging to note that astrophysical data is rapidly reaching the point 
where it can put severe constraints on large classes of inflationary models. In particular, 
the recent analysis by Seljak et al. S uses a large set of astrophysical data sources to 
argue for tight constraints on the primordial spectrum, which rule out many “standard” 
models of inflation. In particular, Ad 4 is excluded at the 3cr level, while the single-field 
m 2 (f) 2 model makes predictions that differ from the central value at the 2 a level. Obviously 
these constraints on the primordial spectrum will need to be confirmed (and will hopefully 
tighten) as further datasets become available. As noted earlier, we have not fully analyzed 
the inflationary predictions of N-flation in the general (3^0 case, and we do not compute 
direct astrophysical constraints on the parameter values. However, the overall size of the 
error ellipses in [52[] suggests that observations may soon put tight constraints on the N- 
flation parameter space (if they do not do so already), providing further motivation for a 
careful analysis of the cosmological predictions of this model. 


9. Possible Exceptions to our Results 

In this section we collect the conditions and assumptions that were required in order 
to ensure the controllability of N-flation generally and to derive the Marcenko-Pastur 
spectrum. 

First, in §9.1, we review the constraints implied by controllable perturbative and 
instanton expansions. These requirements, which were stated in (1J, have little to do with 
the inflationary dynamics per se; instead, they should be thought of as necessary conditions 
for a controllable trans-Planckian vev created by N string axions. 

Then, in §9.2, we summarize the assumptions and constraints that apply particularly 
to our random matrix ensemble of N-flation models, with the intention of revealing the 
exceptional classes of N-flation models that might avoid our conclusions. 
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9.1. Conditions for a Trans-Planckian Displacement 

One of the main results of [;1[] is the construction of plausibly trans-Planckian vevs 
A$ > M p from multiple string axions. As explained in detail in |jT| . radiative corrections 
diminish the effective size of A$ and reduce the expected scaling 


(A$) 2 cx N 


(9.1) 


to 


(A$) 2 oc 


N 


X(M) 


(9.2) 


where we recall that in our notation, P — 1 is the number of complex structure moduli, 
and 

x{M) = 2iV — 2(P — 1) (9.3) 

The result is that A$/M p can be large, but only when there are partial cancellations 
between N and P. This is a powerful reason to prefer models with (3 ~ where this 
cancellation is almost total. (An unfortunate consequence of this numerical, rather than 
parametric, success is that factors such as 27r are of considerable importance in determining 
the output of a given model.) 

Furthermore, although one might hope to construct trans-Planckian vevs in a string 
vacuum with extremely weak coupling and extremely large volume, this has not yet been 
achieved. The problem is that the axion decay constants diminish as the volume increases, 
so it is not possible to work at arbitrarily large volume. 

Let us therefore examine the conditions for control of the two perturbative expansions 
of string theory: the string loop expansion, controlled by g s , and the sigma model expan¬ 
sion, controlled by a.'/R, where A is a typical length-scale in the compact space. Except 
when noted, we are here simply summarizing the arguments of [|T]. 20 

The leading contribution to the Planck mass from sigma model corrections is 


dg'Mp x (M)C(3K 3 
8tt 3 f/ 6 


M 2 

- LV - L p 


(9.4) 


which depends on N through y(M). Next, the dangerous corrections to M p from the string 
loop expansion are those that have the leading N-scaling, namely 


8g.Ml 

Ml 


V 167T 2 7q 


M, 


kk 


M n 


2 q 


(9.5) 


20 


We thank J. Maldacena for very instructive discussions of these arguments. 
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In the final expression, we have included the usual field theory loop factor of 167 t 2 , as well 
as a parameter that measures the deviation of the actual prefactor, in the string theory 
q-loop result, from this naive expectation. (In [1], the very conservative choice ^ 

was made. It would be worthwhile to compute this factor directly for the one-loop case.) 
The ratio of the Kaluza-Klein mass M^k to the Planck mass arises in ( |9.5|) for the following 
reason. We are considering loop corrections to the four-dimensional Planck mass that arise 
from axions, or their moduli partners, circulating in loops. These fields have couplings of 
gravitational strength, and the diagrams must be cut off at the compactification scale; 
combining these results leads to ( |9.5| ).il 

The minimal requirements for perturbative control are 


5 a' Ml 
M 2 

- LV - L p 


< i 


(9.6) 


Ml 


< 1 


(9.7) 


which imply constraints on g s and on the compactification volume. Among these is the 
important result that ( |9.6|) is most readily satisfied when (3 ~ 

In addition, we must ensure that the single-instanton contribution to the superpoten¬ 
tial is a good approximation. This follows from 


Pi > 1 


i = 1...N 


(9.8) 


A skeptical reader may be concerned that we have dropped multi-instanton terms in the su¬ 
perpotential ( |2.5|) but retained the products of two single-instanton terms in the potential. 
This is not inconsistent: the leading contributions to the potential from a multi-instanton 
term in the superpotential have the schematic forms 

AF oc D a W 0 e~ 2n(Ti+Tj) + c.c. (9.9) 


where a is a complex structure modulus, and 

Ay oc W 0 e~ 2n ( Ti+T ^ + c.c. (9.10) 


21 In [jlj, higher-loop terms were conservatively taken to be proportional to a single power of 

M 2 

-jjT- Even in this case, the requirements of control can be met. 

ivtfp 
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The factor D a Wo in the former would vanish |[4[| if the background fluxes were imaginary 
self-dual. In a vacuum stabilized by nonperturbative effects, this factor will be nonvanish¬ 
ing, but nonperturbatively small J53|. Thus, (|TDj) is negligible compared to the terms we 
have retained. Next, ( |9.10| ) is highly suppressed because of the smallness of Wq: to set an 
appropriate overall mass scale, we needed to work with a flux configuration in which Wq is 
very small [|37|], exactly as in Finally, notice that our random matrix arguments could 
be extended even to cases in which multi-instantons are relevant. 

It was argued in |JJ that the intersection of the constraints mum (Ei) is nonempty 
and includes configurations with trans-Planckian displacements. Even so, verifying this by 
direct computation in an explicit example would be most valuable. 


9.2. Conditions for an N-flation Model with Marcenko-Pastur Spectrum 

For definiteness, we have focused our attention on KKLT vacua in type IIB string the¬ 
ory. However, scenarios other than that of KKLT have some notable virtues for N-flation. 
Type IIA moduli stabilization with fluxes |fTT| produces a hierarchy between the scale of 
moduli stabilization and the scale at which nonperturbative effects produce axion masses. 
This reduces the challenge of arranging that the Kahler moduli f); are massive enough to 
remain at their minima during inflation. Heterotic string vacua have the advantage that 
the constraints of perturbative control reviewed in §9.1 can be less stringent for NS axions 
than for RR axions, because the relation between the compactification volume and the 
decay constants is different, by factors of the string coupling, in the two cases. The axions 
suitable for N-flation in known, stable type II vacua are RR axions whose nonperturba¬ 
tive potentials come from D-brane instantons, whereas in the heterotic string the relevant 
axions come from the NS B-field and receive their potentials from worldsheet instantons. 

More importantly, for any vacuum in the above classes, the only essential property 
that we provided as an input to the supergravity potential ( |T2| ) is that the axion potential 
arises from F-terms that vanish when the axions are at their origin. As explained in §2.3, 
in practice we require only the natural assumption that any F-term energy present at the 
axions’ minimum is negligibly small compared to the energy scale of inflation. In this 
sense, our results are extremely general: we expect that the Marcenko-Pastur spectrum of 
masses is characteristic of N-axion potentials in most of the approximately F-flat vacua in 
any supergravity theory.Hl 


22 Although our primary motivation was inflation, 
potential, whether or not this is relevant for inflation. 


I) is truly a result about the N-axion 
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Axion inflation around a vacuum that is not even approximately F-flat could lead to 


a mass matrix whose form is different from (|4.7|) . However, we expect that in such a case 
the mass spectrum will be an appropriate mixture of the Wigner semicircle and Marcenko- 
Pastur distributions. Given a specific, well-motivated model that is not F-flat, one could 
easily apply our methods to determine n s . 

We approximated the entries (|4.2| ) as independent, random variables. Strong corre¬ 
lations among these terms would violate this assumption. However, it was shown in [|54 


that the Wigner law persists in ensembles with substantial correlations among the matrix 
entries, and it is reasonable to anticipate a similar result for the Marcenko-Pastur law. We 
are unaware of any motivation for a correlation that is substantial enough to invalidate 
( |4.8| ) . Of course, it would be extremely interesting to find exceptions; the effect of known 
correlations could possibly be analyzed through methods similar to those of |)5|. Moreover, 
even in the presence of strong correlations and totally non-random behavior, some basic 
results are available. Using only the fact that M ( fS.fJ| ) is positive-definite, and making no 
assumptions whatsoever about the statistical properties of its entries, it follows |5(J] that 
for any i ^ j, 


m„ 


rn„ 


> 2|Af, 


ij i 


(9.11) 


where the l.li.s., the difference between the largest and smallest eigenvalues of M, is known 
as the spread of M. This implies a weak but unavoidable lower bound on the difference 
between the largest and smallest axion masses-squared. 

We did not have to assume that the distributions are Gaussian or otherwise simple, 
or that the mean /i is zero, because the result of Marcenko and Pastur can be shown m 
to follow without these conditions. We also did not have to take the N —> oo limit. The 
Marcenko-Pastur law applies as N —> oo, but is an excellent approximation for the values 
of N > 200 that were already required for N-flation to be successful. Corrections to the 
moments ( |6.14|) are of order -h [flip and are too small to be relevant in our results. 

We did have to make a restrictive assumption (as was done in [IX| for simplicity, 
though some statements made there do apply to general initial conditions) about the 
initial conditions, namely that the initial vevs obey fa < Cl Vf. The motion of N axions 
with large vevs is governed by ( ft.l(J| ), but we have left a study of the dynamics of such a 
system as an interesting problem for the future. 

We noted that, as a result of the renormalization of the Newton constant explored in 
[I[] and reviewed in §9.1, models with (3 ~ \ are preferred. In such models the spectrum 
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is typically more red, for a range of initial conditions, than in single-field m 2 (j) 2 inflation. 
An obvious way to avoid this constraint would be to demonstrate that models with much 
smaller values of /? are in fact controllable and can give enough inflation. 


10. Conclusions 


We began with the observation that substantial cross-couplings between axions are a 
generic feature of N-axion models in string theory. Moreover, complicated kinetic terms 
further connect the various axions. Finally, as suggested in [dj, the masses of the individual 
fields are not identical, in general, but are spread over some range. Each of these effects 
can lead to nontrivial departures from the simplest realization of assisted inflation: the 
fields do not roll independently, and the dynamics is more complicated than that of chaotic 
inflation with the potential m 2 (j) 2 ||. 

When the initial axion vevs are taken to be small in units of the appropriate decay 
constants, the cosine potentials may be expanded, and the leading term in the resulting po¬ 
tential is quadratic. In this setting, one may rotate and rescale fields to produce canonical 
kinetic terms and a diagonal mass matrix. The result is a model specified by a spectrum 
of masses for N uncoupled axions. 

In this basis in which the axion kinetic terms are canonical and the mass matrix is 
diagonal, the dynamics appears relatively simple, though not quite as simple as in ordinary 
m 2 (j) 2 inflation. The most massive axions fall to their minima first, diminishing the Hubble 
friction and releasing progressively lighter axions to fall to their minima. The Hubble 
parameter, and hence the spectral index of the scalar perturbations, therefore depends on 
the spacing in mass between adjacent axions. In any given compactification, this spacing 
depends on a formidable array of microscopic details. However, the statistics of the spread 
in masses can be characterized much more simply. 

We used random matrix techniques to cornputell the typical axion mass spectrum, 
taking as input the structure of the supergravity potential. We found that the spectrum 
is given by a simple and universal analytic formula ( |4.8|) , the Marcenko-Pastur law ||39|| , 
in the limit N —> oo. For realistic values of N. the agreement between this prediction and 
the results of Monte Carlo simulations of the mass spectrum is excellent. 


23 The overall mass scale m must be fixed by requiring the correct amplitude of density pertur¬ 
bations, as in jl|. 
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The shape of the Marcenko-Pastur curve (~) depends on a single rational parameter 
ft ( |4. 11 1 ) . which is determined by the dimensions of the Kahler and complex structure 
moduli spaces. Our results and constraints are therefore phrased in terms of ft alone. This 
is a rather dramatic simplification from the naive expectation that a randomly chosen 
model of the form ( p.3|) would be characterized by O(N) independent parameters. Even 
worse, in the full model ( |3.10| ) in which the initial displacements are so large that terms 
beyond quadratic in the axions must be kept, at least 0(N 2 ) parameters would be required. 
In this sense, random matrix models of N-flation are surprisingly predictive. 

We then studied the inflationary trajectory of a general N-flation model. Using the 
distribution Q4.8Q to generate an ensemble of mass matrices simulating an ensemble of N- 
flation models, we numerically integrated the equations of motion for this system. Next, 
we studied the adiabatic perturbations in a general N-flation model. We gave analytical 
expressions for the scalar power spectrum ( |7.4| ) and its tilt ( [7.10| ) and running ( |7. 1 from 
which, given any spectrum of axion masses and choice of initial conditions, one can readily 
deduce the cosmological observables. 

We found that the scalar power spectrum resulting from a generic N-flation model 
is distinguishable, in principle, from the single-field m 2 (f) 2 inflation result. In the case of 
the initial conditions that we argued are most plausible, the spectral index can be written 
very simply in terms of 3. as ( |7.17| ) . For nonzero ft- the spectrum is more red than in the 
associated single-field model; for the preferred value ft « \ the change in n s — I is around 
25%. Let us stress here that there do exist parameter regimes, e.g. ft —> 0, in which an 
N-flation model is observationally similar to m 2 (ft 1 inflation. However, to the extent that 
these are unusual parameter choices, we find that N-flation ‘usually’ differs from m 2 (j) 2 
inflation. 

Because the dynamics is not independent of the initial conditions, even at late times, 
we presented our general results as functions of the initial conditions. For example, knowing 
only ft and not the initial conditions, we cannot predict the exact tilt of the scalar power 
spectrum, but we can make the general prediction that the spectrum will be more red than 
in the single-field case, by an amount that depends on ft. In reasonable special cases in 
which the field vevs or the energy content of each field are roughly equal, the results are 
substantially simpler. In this connection, it would be interesting to characterize the most 
plausible initial conditions for models with large numbers of moduli. 

Although we propose that random matrix theory is a very useful tool for analyzing N- 
flation, we certainly do not claim that random matrix results produce or reveal a flattening 
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of the inflaton potential. These techniques merely allow us to extract a nearly-universal 
axion mass spectrum in a context where direct computation of the masses is unmanageable. 
In fact, our results imply that the scalar power spectrum of a randomly chosen N-flation 
model is typically more red than in m 2 (f) 2 inflation. Thus, we actually find that a generic 
N-flation potential is fractionally steeper than in the simplest model proposed in [Q. For¬ 
tunately, this slight steepening does not spoil slow-roll inflation, and, as we explained in 
§8, it can sometimes yield a testable prediction. 

One interpretation of our results is that, in certain circumstances, assisted inflation 
has a novel advantage over the a priori much simpler single-field inflation models. Suppose 
one chooses an embedding of inflation in string theory that involves one or a few dynamical 
fields. The exact form of the inflaton potential will, of course, depend on the specifics of the 
compactification. Given an ensemble of realizations of this model, with differing choices 
of the microscopic parameters, the mean behavior will give only a qualitative guide to 
the behavior of a particular realization: fluctuations from the mean are not suppressed by 
large-number effects. Consider, in contrast, an ensemble of realizations of N-flation, or of 
some other sort of assisted inflation. The large-N simplifications of random matrix theory, 
as well as the well-known suppression of fluctuations from the mean, ensure that a typical 
realization of N-flation is very well described by the mean behavior in this ensemble of 
realizations. This greatly enhances the predictivity of this class of models. In a sense, this 
runs counter to the naive argument that the complexity of the string landscape makes it 
almost impossible to draw conclusions about the inflationary models that exist within it. 
One might seek other situations where random matrix theory can be combined with the 
large dimensionality of the landscape to draw generic conclusions about classes of models 

MM 

In particular, it would be very interesting to find other examples in which assisted 
inflation becomes somewhat more predictive and less dependent on microscopic data as 
a consequence of large-N statistical effects. Not every assisted inflation model is usefully 
specified by a matrix or a set of matrices, so random matrix results per se are not always 
applicable. Even so, for a sufficiently large number of fields, statistical reasoning may still 
produce useful simplifications. 

We focused, for definiteness, on KKLT compactifications of type IIB string theory, 
but as explained in §9, it is entirely clear that our conclusions generalize to other contexts. 
Even so, it would be valuable to examine the detailed properties of the N-flation potential 
in other string theories or in M-theory. 
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We have left untouched many other interesting questions about N-flation. In particu¬ 
lar, although we have learned something about the statistical properties of an ensemble of 
N-flation models, we have not provided anything like a concrete realization of the model, 
e.g. in a specific orientifold. This would be a worthwhile endeavor. 
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Appendix A. Statistical Suppression of Cross-Couplings 


In this Appendix we will compare the sizes of diagonal and off-diagonal elements in the 
mass matrix ( |3.9|) . We will show that the diagonal elements M it are normally distributed 
around a mean that is larger, by a factor ~ a /N + P, than the expected magnitude (|\) 
of off-diagonal entries My, i ^ j. 


As we explained in §4, the basic structure of (|3.9|) can be captured by 


Mij — R,aR 


(A.l) 


where R is a random (N + P) x N matrix that is related by a rotation to Haj = DaCj. 
The entries of R may have either sign. As we demonstrated in §4.1, the entries of R have 
a mean // that is small compared to the standard deviation a. In light of this property, we 
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assume 


a 2 . 


for simplicity® that the R r/ \ are normally distributed around zero with variance 

The expected value of a diagonal element Mu is immediately seen to be 

{Mu) = (N + P)a 2 (A.2) 

The off-diagonal elements are only slightly more involved. We start by examining 


M^ij ^ ^ RiARjA, 


(A.3) 


which contains a sum of N+P products of distinct elements Rj ,\ x R-ja- Using the fact that 
the distribution of a product of normal distributions, each with variance a 2 , has variance 
a 4 , and applying the central limit theorem, we find that the My are normally distributed 
with zero mean and variance (N + _P)cr 4 .i 1 Thus, we find 


(IMA) = ^AnTpc 


(A.4) 


It follows that the ratio of the expected magnitude of the diagonal terms to that of the 
off-diagonal terms is 


( 


Mi, 


I Mi 


-} = y/WTP 


ij l 


(A.5) 


which is, up to a factor very near unity, the scaling we found in our Monte Carlo simulation. 
Thus, given the knowledge of the form of matrix product dictated by the supergravity 
potential, and using our result from §4.1 that the D aC[ have mean much smaller than one, 
we have derived the bimodal distribution found in §5, including the suppression of 

off-diagonal elements. 


Appendix B. Inflation in the Presence of Cross-Couplings 

In this Appendix we will study the inflationary dynamics of N canonically-normalized 
fields coupled by a general quadratic potential, in some basis in which the mass matrix is 
not diagonal. 

24 It is straightforward to extend our arguments, though not necessarily the order-one coeffi¬ 
cients, to much more general distributions. 

2o Because we consider a sum of N terms, not a mean, the scaling is superficially different from 
the classic application of the central limit theorem. 
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Let us parametrize the mass matrix as 


Mij — m 2 5ij + em 2 Ejj (B.l) 

where £ is a dimensionless constant, m sets the average mass scale, and E,; ? - is a symmetric 
matrix whose entries have mean zero and variance a 2 . 

We would like to determine the critical value of £ for which the dynamics is dominated 
by the diagonal terms, with no substantial corrections due to cross-couplings. Let us first 
compute the total potential V when each of the fa is near some average value &. We find 

V ~ (n ± eNaj (B.2) 

where the factor of N in the final term comes from the random walk of 0{N 2 ) entries in 
E. (By ±£iVa we are indicating the typical range of variation of V.) 

We now observe that 


2 diV — 2m 2 fa + £m 2 'Eij(j)j ~ 2m-^ ± £m 2 aVN(f> 


(B.3) 


where we have noted that the final term is a sum of N terms distributed around zero with 
variance a 2 . 

With these results we can read off the slow roll parameters: 


e, = 


°i vi 


Vii = 


Vij = 


m 2 (j) 2 N 


2 m 2 M 2 
fh 2 (j) 2 N 


e^jM 2 

4> 2 N 


ay/N<f>\ 2 n n,—2 

r ) (1 ±£(t) 

(B-4) 

(1 ± £cr) _1 

(B.5) 

(1 ± £cr) -1 

(B.6) 


Clearly £cr 1 suffices to arrange that ( |B. 5|) is unmodified by the presence of the 
cross-terms E. Moreover, because (|E^|) ~ a, this same condition on £cr implies that 
is negligible compared to rja. The most stringent constraint is that the cross-terms do not 
substantially change e. This requires £cr < W -1 / 2 . 

We conclude that £cr < iV -1 / 2 is the necessary condition for the slow-roll parameters 
to be unmodified by the presence of quadratic cross-couplings E,;,-. However, we just 


showed above that in the case of N-flation, eo 


VN+P- 


Thus, the cross-couplings in the 
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quadratic-potential limit of N-flation are naturally small enough to leave the inflationary 
dynamics unaffected. (We expect that cross-couplings cannot be neglected in the full 
N-axion potential, as explained in §3.2.) 

Two comments on this issue are in order. First, the fact that the cross-terms are 
dynamically unimportant does not mean that they should simply be omitted: the dis¬ 
tribution of eigenvalues of the mass matrix is different from the distribution of diagonal 
elements. Second, we interpret the suppression of cross-couplings as originating in statis¬ 
tical considerations, not more elementary physical ones, as evidenced by the fact that the 
suppression factor £ is the square root of an integer, and not a ratio of mass scales. 
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